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ABSTRACT 



We used fully cosmological, high resolution N-body + SPH simulations 
to follow the formation of disk galaxies with rotational velocities between 135 
and 270 km/sec in a ACDM universe. The simulations include gas cooling, 
star formation, the effects of a uniform UV background and a physically 
motivated description of feedback from supernovae. The host dark matter 
halos have a spin and last major merger redshift typical of galaxy sized halos 
as measured in recent large scale N-Body simulations. The simulated 
galaxies form rotationally supported disks with realistic exponential scale 
lengths and fall on both the I-band and baryonic Tully Fisher relations. An 
extended stellar disk forms inside the Milky Way sized halo immediately 
after the last major merger. The combination of UV background and SN 
feedback drastically reduces the number of visible satellites orbiting inside a 
Milky Way sized halo, bringing it in fair agreement with observations. Our 
simulations predict that the average age of a primary galaxy's stellar 
population decreases with mass, because feedback delays star formation in 
less massive galaxies. Galaxies have stellar masses and current star formation 
rates as a function of total mass that are in good agreement with 
observational data. We discuss how both high mass and force resolution and 
a realistic description of star formation and feedback are important 
ingredients to match the observed properties of galaxies. 

Key words: galaxies: evolution — galaxies: formation — methods: N-Body 
simulations 



1 INTRODUCTION 

In a universe dominated by Cold Dark Matter ( CDM) 
and a cosmological const an t ([Blumcntl ial et al.l 



19841: iBardeen et al.l 1 19861 : ISpergel et alj l2003l : 



Peebles fc Ratral |2003| ). galaxy formation and evo- 
lution is a complex combination of hierarchical 
clustering, gas dissipation, merging and secular evolu- 
tion. While gravit y drives the botto m-up assembly of 
cosmic structures (|Davis et al.lll98^ . gas cools at the 
centers of dark matter halos and acquires angular mo- 
mentum through tidal torques from nearby structures 



(|Fall fc Efstathioulll980l : lFaiilll983D . Cool gas eventually 
fragments due to Jeans instabiliti es and forms stars 
(|White fc R"e5ll978l : [Peebles! Il969'l. Gravitational in- 
teractions between galaxies can insti gate star formation 
and transform disks into spheroids I Toomre fc Toomrd 
Il972l ) and may dominate galaxy evoluti on at early 
epoch s especially in dense environments l|Frenk et al.l 
Il985l ). Secular processes, such as external gas accretion 
or gas displacement due to bars, may play a more 
i mportant role in ga l axy e v olution at later e pochs 
jValenzuela fc Klypin' '2003'; 'Pebattis ta et al.l |2004| : 
iKormendv fc Kennicutt .2004 1 . Galaxy formation in 
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this framework is a continuous, ongoing process where 
the observable properties of galaxies are a function 
of t heir merging his t ories, masses and environments 
(e.g. lShen et a l 2003': ' BlantonI [20031 : iBrinchmann et all 
l2004l : lBund~ ;t al.„2005f r 

N-Body/gasdynamical simulations have become 
the primary tools with which to model galaxy forma- 
tion in a cosmological context. They are necessary to 
follow the evolution of the internal structure of galaxies 
as well as the compl ex interplay between b aryon 
cooling and 
19921: 



Ha 



19961: fO uinn ct al. 
Mac Low & Fcrraraj 



feedback llLake fc Carlberd Il988l: iKatz] 
TO fc White^ Il994l; [Barnes fc Hernquist 



Thacker fc CouchmanI 



Haehnc lt et al. 1998 ; 

Cen fc Ostrikerl Il999l: 

2000l: iNavarro fc Steinmet j 



20001: iTittlev et all l200ll : ISommer-Larsen et al.l l2003l : 

Governato et al.l2004 '). However, numerical simulations 
have had difficulties reproducing the properties of 
real galaxies. Early work reported a catastrophic loss 
of angular momentum in the baryonic component of 
simulated galaxies, leading to the formation of galaxies 
dominated by a central concentration of cold baryons 
l|Navarro fc White 1994 ). As p r edicted by early the- 
oretic al models l| White fc Reed Il978l : IWhite fc Freii^ 
I1991I ). simulations without strong energy feedback 
from stellar processes have also suffered from overly 
efflcient star formation. This causes galaxies to form 
most of their stellar component as soon as gas is able 
to collapse and cool within small halos , well before 
the assembly of the main progenitor (jBalogh et al.l 
I2OOII ). Dense clumps of stars and cold gas might then 
spiral rapidly at the ce nter of the main galaxy du e 
to dynamical friction ([Navarro fc Steinmetj 1 19971 ). 
These shortcomings resulted in compact galaxies with a 
central excess of old stars^and offset from the o bserved 
TuUy-Fisher relation (jGiovancUi ct~al] Il997l '). The 
offset was ^ 1.5 mag for an = 1 CDM cosmology 
([Navarro fc Steinmetzll2000l: [Eke et al.ll200l[ ') and ~ 0.5 



mag in a high as ACDM ([Eke et al.ll20oi[ r 

More recently, attention has been drawn to the 
"missing galaxy problem". Dark matter only simula- 
tions predict over an order of magnitude more subha- 
los around Milky Way-like galaxies than the number 
of dwarfs observed around the Milky Way and M31 
jKlvpin et al. [ 1 19991 : [Moore et al.l Il999l : IWillman et al.l 
[20041 ). Observations also show that more massive 
late-type galaxies have older s tellar populations than 
less massive late-type galaxies (iMacArthur et al] [20041 : 
[Ferreras etakl [20041 : [Callazzi et al.[ [20051 ). in possible 
contradiction with the fact that in CDM scenarios less 
massive halos assemble (on average) first. 

The angular momentum and missing galaxy prob- 
lems have often been linked to two major stumbling 
blocks for numerical galaxy formation: i) insufficient 
numerical resolution and ii) inaccurate treatment of 
feedback due to star formation. Only the most recent 



gas dynamical simulations have achieved spatial resolu- 
tion sufficient to resolv e th e disk scale lengths of disks . 
[Governato et al.[ ([2004[ ) and [Kaufmann fc Mavel ([2006[ ) 
have shown that poor mass and spatial resolution might 
lead to significant numerical angular momentum loss in 
baryonic disks embedded in dark matter halos. These 
works suggested that at least 10^ DM particles (and 
equivalent numbers in gas and star particles) within the 
virial radius as well as force resolution of < 1 kpc are re- 
quired to faithfully simulate the formation and dynam- 
ical evolution of disk systems over cosmic times. Re- 
cent simulations ("Brook et al.' '2004'; 'Abad i et al.[ [20031 : 
Robertson et al. 2004; Okamoto ct al. 200i) that used 
a large number of particles and a high force resolution 
have indeed formed galaxies with extended stellar disks. 
However, most of the galaxies formed in most of these 
simulations also had a massive spheroidal component or 
disks only partially supported by rotation. 

The impossibility of directly resolving the scales at 
which star formation and feedback happen (a few pc) 
makes it necessary to develop simplified models to de- 
scribe star formation and subsequent energy feedback 
from supcrnovae (SNe) at gala ctic scales (0.1 - 1 kpc) 
(Efstathiou 2000; Silk 200 ll: [Wada fc NorinanI I2OO1I: 
Ferrara ,2002 : .Krumholz fc McKed [2005[ : [Slvz et al.[ 
2005h . In early simulations, thermal energy that was 



simply added to gas surrounding star forming regions 
was q uickly radiated away, resulting in overcooled gas 



([Katj[l99j ). However, heated, diffuse baryons result- 
ing from a proper treatment of feedback may be less 
susceptibl e to angular momentum t ransfer to halo 
particles ( Mo fc Miralda-Escudc 199^; [Eke et al.[ [2OO0I : 
[Mailer fc DekeL,2002: ,p'Onghia.,200Q). To address this 



Yepes et all ([1993) and 



shortcoming. iGerritsenI (Il997l ) , ., . . 

[Thacker fc CouchmanI ([2OO0I ) modeled the pressure sup- 
port stimulated by gas turbulence by shutting off cool- 
ing on timescales of a few million years, and r eported 
the formation of larger di sks. Brook et al. used a 

scheme similar to that of [Thacker fc CouchmaiJ (|2000l ) 
and also found it produced reasonably realistic galaxy 
disks when applied to simplified cosmological initial 
conditions. [Robertson et al.l ([2004[ ) used an approach 
similar to that of [Springell (|2000l ) and found that the 
pressure support produced by their multiphase treat- 
ment of the ISM was an important factor in the for- 
mation of large disks. Supernova feedback and the UV 
background, being able to reduce gas retention and 
accretion in halos with low virial temperature, may 
also s olve the problem o f the overabun d ance of satel- 
lites (O uinn et all il996': 'Cnodin 2000; Bens on et al 
|2002. : Dekel fc Wocli2003 : .Monaca.2004 : Kravtsov et al 
I2OOJ) . Other rece nt papers have iriyestig ated the roles of 
a top heavy IMF ([Okamoto et al.|[2005l ) a nd of massive 
galaxy outfiows driven by QSO activity dBinnev et al.[ 
[200ll : lGranato et al.|[20o3 : [Di Matteo et al.[[2005^ in reg- 
ulating galaxy formation. 
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In this paper, we study the effect of feedback on the 
structure and satellites of disk galaxies formed within 
cosmological halos spanning a significant mass range 
(from 10^^ Mq down to 10*). We improve over previous 
work in two ways. First, the mass and spatial resolution 
of these simulations are sufficient to resolve (i) the struc- 
ture of present day disks without being significantly lim- 
ited by resolution and (ii) the sub-halo population for 
each galaxy in our sample down to circular velocities of 
about 20% of their parent halo allowing us study the 
basic properties of galactic satellites. Second, we use a 
revised implementation of the star formation algorithm 
introduced byJKatz (19 9j) combined with the feedback 
rec ipe introduced bvlGerritsenI (Il997l ) and then studied 
by iThacker fc CouchmanI ||2000|~ Details of our imple- 
mentation can be found in Stinson et al.l l|2006l ) (here- 
after S06). In this paper we calibrate the free parameters 
of our algorithm to accurately describe the star forma- 
tion in isolated galaxy models of present day galaxies. 
We then apply this algorithm to cosmological simula- 
tions of individual high resolution galaxies spanning a 
decade in mass and analyze their properties. 

We focus our analysis in this paper on three funda- 
mental properties of present day galaxies: 

• The abundance and luminosities of galactic satel- 
lites. 

• T he Tully-Fisher and baryonic TuUy-F isher rela- 
tions (|GiovaneUi et al.lll997l : lMcGaughll2005l ). 

• Global star formation histories (SFHs) and z = 
star formation rates. 

In §2, we briefly describe the code and the star for- 
mation algorithm. A detailed description of our star for- 
mation and feedback algorithm with additional tests to 
those presented here is presented in S06. In §3, we sum- 
marize the cosmological runs performed with a range 
of feedback algorithms, and in §4-9 we discuss and de- 
scribe in detail the properties of our simulated galaxies. 
We plan to explore other aspects of both low and high 
redshift galaxy formation in future papers. 

2 THE CODE AND THE STAR FORMA- 
TION ALGORITHM 

2.1 GASOLINE 



We used GASOLINE l|Wadslev et al.ll2004l ). a smooth 
particle hydrodynamic (SPH), parallel treecode that is 
both spatially and temporally adaptive with individual 
time steps. The version of GASOLINE we used imple- 
mented i) Compton and radiative cooling, ii) star for- 
mation an d a supernova feed back as described in detail 
in S06 and lKatz et al.l (| 19961 ). and iii) a UV background 
following an updated version of lHaardt fc Madau l|l996l ) 
(Haardt 2005, private communication) starting at z = 6. 
The opening angle, 9, was 0.55 until z = 2 and 0.75 



there after, and the time-step ping criterion, i], was 0.2, 
as m iDiemand et~aLl (|2004( '). The minimum smooth- 
ing length allowed is equal to 0.1 times the gravi- 
tational softening. Exactly 32 neighbors are used for 
SPH calculations. An ideal gas of primordial compo- 
sition is assumed. We calculate the cooling/heating 
rate and ionization state of each particle by assum- 
ing coUisional ionization equilibrium and the presence 
of the time dependent, but uniform UV background. 
The code tre ats artificial viscosity as suggested in 
lBalsaral(|l997l). The energy equation is c o mputed asym- 
metrically fepr ingel fc Hernguisd |2002| : |e vrar 

3 1 19881 : 

jVIonaghan 1992). This approach avoids the energy con- 
servation and negative energy problems of the arith- 
metic and geometric implementations and converges to 
the high resolut ion answer as well as other recently pro- 
posed meth ods ('Bcnj |l990l : ISpringel fc Hernquistli20o3 : 
I Wadsley et al. 2004). 

2.2 The Star Formation and SN Feedback Algorithm 

Our star formation (SF) algorithm can be broken down 
into three main parts: (i) identifying the star forming 
regions, (ii) forming stars and (iii) treating stellar evo- 
lution including such effects as mass loss, SN winds and 
metal enrichment. The algorithm has only three free pa- 
rameters: c*, eSN, and (3. We define these parameters 
and summarize the main features of the algorithm be- 
low. This algorithm is described in detail in S06. 

2.2.1 Criteria for Star Formation 

The criteria for a gas particle to become eligible for 
star formation are: (i) Its temperature is colder than 
30,000 K; (ii) Local gas density is higher than 0.1 cm~^; 
(iii) The gas particle is part of a converging flow mea- 
sured over the 32 nearest neighbors. We do not re- 
quire the criterion th at gas particles are Jeans unstable 
(|Kawata et al.ll2005l '). A minimum gas overdensity Sp/p 
is also required to avoid spurious star formation at very 
high z. We base the rate at which gas is converted into 
stars, dpi,/dt, on the relation 



dp* 
Itt 



^ Pa 



3/2 



(1) 



where p represents the volume density. Using the fact 
that dynamical time, tdyn oc p~^^^, we express this as 

dp* Pga 
dt td 



(2) 



yn 



where we have introduced a constant efficiency factor c* 
to enable us to calibrate the star formation algorithm 
to match star formation rates observed for z = galax- 
ies. The mass of star particles formed is fixed to 30% 
of its parent gas particle initial mass (|Okamoto et al.l 
120051 '). Once the particle passes the above criteria, to 
implement Equation (1) in a discrete system we assign 
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a probability p that a star will actually be spawned from 
its parent gas particle: 

MsF 



-(1 



(3) 



Mgp 

where we have introduced Msf, the spawning mass for 
star particles, Mgp, the mass of the gas particle that is 
creating the star, At, the star formation timescale (one 
Myr in all of the simulations described in this paper) 
and tform, which is either the dynamical time or the 
cooling time, whichever is longer. Gas particles in dense 
regions with shorter dynamical times will form stars at 
a higher probability. 

2.2.2 Feedback and Metal Enrichment From Super- 
novae 

Our implementation of feedbac k qualitatively follows 
the a lgorithm implemented by Ixhacker fc CouchmanI 
(|2000l ). We assume that the energy released into the ISM 
turns into turbulent motions (at unresolved scales) and 
is partially dissipated, preventing the gas from cooling 
and forming stars. We determine the number of Super- 
nova type la fc II th at occur during each time step from 
iRaiteri et al] (119961 ) using a Miller-Scalo IMF and the 
stellar lifetimes of stars. We then multiply the number of 
both SN types that explode by a fraction of the canoni- 
cal 10*^ ergs/SN times a fixed efficiency term (eSN), and 
distribute that energy to the surrounding gas particles. 
In §2.3.2, we explore efficiency values in the range 0.1 
to 0.6. At an efficiency eSN = 0.1 and with the adopted 
IMF, 7.65 X 10^*^ ergs of energy are deposited into the 
surrounding gas for every one Mq of star formed. En- 
ergy is distributed using the smoothing kernel over the 
32 nearest neighbor particles. In our algorithm the time 
scale for the cooling shutoff and the amount of mass 
affected are physically motivated and chosen follo wing 
typical values from the iMcKee fc Ostrikeij (Il977l ) and 
lOstriker fc McKed (|l988t ) blast wave model. When en- 
ergy injection comes from a SNII event we disable the 
radiative cooling for 30 million years in a number of the 
nearest neighbor particles that satisfy: 



PMsNII > 



~ pave 



(4) 



where Msnii is the mass of supernovae produced in a 
star in a given timestep, r is the distance from the star 
to the gas particle in question and pave is the local gas 
density and /9 is a normalization factor. In §2.3.2, we 
explore /3 values in the range 0.05 to 0.2. For each SF 
event the maximum number of particles that can have 
the cooling disabled as from eq.(4) is the number of SPH 
neighbors (32 in our simulations). 

Once formed, we treat each star particle as a single 
stellar population with uniform metallicity. The code 
also keeps track of mass loss from stellar winds. With 
the adopted IMF, star particles lose up to 30-40% of 



their original mass as their underlying stellar population 
ages. This mass gets distributed to nearby gas particles. 
Metals come from both SN e la fc II. Metal enrichment 
follows ^Raiter i et al.l (|l996l ). Like energy, metals are dis- 
tributed using the smoothing kernel over the 32 nearest 
neighbor particles. We plan to study the metallicity of 
the baryonic component of our simulated galaxies and 
to introduce metal lines cooling in future papers. 

2.3 Using Isolated Galaxy Models to Calibrate the SF 
and SN Feedback Algorithms 

We performed a number of simulations of isolated in 
equilibrium galaxy models to study the effect of different 
combinations of the three free star formation parame- 
ters (c*, j3 and eSN) on the properties of galaxies with a 
quiescent star formation and over a range of masses. Pa- 
rameters were varied within a range of plausible values 
suggested by observational constraints. The star forma- 
tion efficiency parameter c* was varied over the range 
[0.01 - 0.4], the mass factor /3 over [0.05-0.4] and the 
SN efiiciency eSN in the range [0.1 - 0.6] (see Table 1). 
This is a subset of the parameter space explored in S06 
where tests were focused on a Milky Way (MW) sized 
galaxy. However, it is important to extend tests to less 
massive galaxies. Hence in this work we followed a com- 
plementary approach, focusing on the properties of a 
much smaller galaxy and then verifying that the best 
parameters so identified would model realistically star 
formation in a Milky Way sized halo. 

2.3.1 The Galaxy Models 

To explore the effects of our feedback recipe we applied 
it to two galaxy models having different peak circular 
velocities: 220 (MWiso) and 70 km/sec (DWFiso) re- 
spectively (where iso stands for "isolated"). Both mod- 
els were built as equ ilibrium config urations using the 
procedure outlined in^ Springell l|2000l ). Rather than try- 
ing to model specific galaxies, we built models consis- 
tent with the trend for smaller galaxies to have smaller 
bulg e components and a larger cold ga s fraction in the 
disk (lMcGaughll2005l : IWest et al.ll2005l ). The models in- 
clude a rotationally supported stellar and gaseous disk, 
a bulge for the MW model only and a dark matter halo 
component extending to the virial radius. Stars start 
with a Toomre parameter Q=2. Disks are built to be 
stable to bar instabilities. 

The dwarf galaxy has no bulge component and gas 
contributes 50% to its total disk mass. 10^ particles are 
in its dark matter halo, and 1.5 x 10** gas and 2 x 
10* star particles are in its disk. The softening length 
is set to about 0.2 times the disk scale length (Rd ~ 
1 kpc). Gas and stars have the same exponential radial 
distribution. Due to its relatively shallow potential, we 
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Run 




/3 


eSN 


SFR 

Me/yr 


^gas 

km/sec 


^star 

km/sec 


R^/Rd Hot/Cold Gas 
volume ratio 


DWFisoOl 


0.05 


0.05 


0.1 


0.20 


12 


12 


0.23 


0.5 


DWFiso02'' 


0.05 


0.05 


0.2 


0.16 


14 


14 


0.31 


0.6 


DWFisoOS" 


0.05 


0.05 


0.4 


0.12 


17 


16 


0.22 


0.6 


DWFiso04" 


0.05 


0.05 


0.6 


0.06 


35 


22 


0.27 


0.9 


DWFisoOS 


0.20 


0.05 


0.2 


0.18 


22 


20 


0.43 


1.0 


DWFisoOe 


0.20 


0.10 


0.2 


0.15 


20 


16 


0.61 


1.5 


DWFisoO? 


0.20 


0.20 


0.2 


0.13 


18 


20 


0.5 


2.0 


DWFisoOS 


0.20 


0.20 


0.6 


0.03 


36 


28 


0.84 


1.6 



Table 1. Summary of the properties of some of the isolated dwarf runs for a given set of c*, fd and eSN parameters. SFR is 
measured over the whole disk after it has set to a quasi steady rate. (Jgas and astar (the azimuthal velocity dispersions) and 
Rz/Rd ^re measured at two stellar disk scale lengths. The volume ratio between the hot gas (T > 10^ K) vs colder gas 
included all gas particles within 2 kpc from the disk plane. ^ Parameters used for the Cosmological Runs 



expect this galaxy to be fairly sensitive to the details of 
the feedback and star formation algorithm. 

Our Milky Way model (MWiso, the sa me as in S06) 
has th e structural parameters outlined in lKlvpin et al.l 
(|2002l ). so providing an excellent fit to the existing data 
of the matter distribution in our own Galaxy. Our Milky 
Way disk is modeled by 4.5 x 10'' star particles and 10* 
gas particles with a softening of 325 pc. This resolution 
is similar to that of our cosmological runs. The stellar 
bulge to disk ratio is 1:7, and gas contributes about 
10% of the total mass in the disk component. Disk scale 
lengths are 3.5 and 7 kpc respectively for stars and gas 
(|Broeils fc Rheelll997l ). 

2.3.2 Comparison Between Models and Observations 

A set of 25 runs allowed us to study in more detail the 
efi'ect of the three free parameters (c*, fH and eSN) in our 
star formation algorithm on the properties of the DW- 
Fiso model. In this sub-section, we highlight the main 
trends of simulated galaxy properties with each of the 
three parameters to motivate our best choice parame- 
ters. 

Table 1 summarizes the results of 8 representative 
DWFiso runs. The parameters used in these 8 runs all 
produce galaxies with star formation rates (SFRs) that 
range between 0.03 and 0.20 M0/yr and have velocity 
dispersions of their cold (T < 30.000 K) gas compo- 
nent that range between 12 and 36 km/sec. These val- 
ues are all consistent with those observed for nearby 
dwarf galaxies, aside from the 2 galaxies (eSN — 0.6) 
with the largest aqas- The SFR measured from the SDSS 
(|Brinchmann et al]|2004l ) for small galaxies with stellar 
masses around 10^ ranges typically between 0.02- 
0.5 Mp/yr and has a m edian value of 0.2-0.3 Mq/jt 
(|Brinchmann et al] |2004| ) . Observed values for the ve- 
locity dispersion of the cold gas component in disk s 
range in the 10-30 km/sec range (|Pizzella et al.ll2004l ). 

Although all of the representative runs produce rea- 



0.2 0.4 0.6 



0.2 0.4 0.6 

eSN 




Figure 1. Properties of the dwarf galaxy model as a 
function of eSN. The efficiency parameter c* and the mass 
factor /3 have both been set to 0.05. Upper left: Star 
formation rates. The two horizontal lines shows the 95% 
contour for galaxies of st ellar mass lO^'^ M0 from the SDSS 
l|Brinchmann et al.ll20oA . Upper right: Disk scale length 
R^ vs. the disk sc ale height R^. The hor i zonta l line is the 
typical value fromlBizvaev fc Mitronoval l|2002l ') and 
lYoachim fc DalcantonI | |2006|) . Lower Left: The hot /cold gas 
volume ratio as a function of cSN. Lower right: Cold gas 
velocity disp ersion. The ho r izonta l value is the upper value 
of cTgaa from lPizzella et~al ] ll2004h . 



sonable SFRs and aga3, only runs with low c* (0.05 in 
the representative sample) produce disks that are thin 
enough to match observations. The observed ratio of 
the stellar disk scale height (Rz) to stellar disk scale 
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Figure 2. The gas projected density for the dwarf galaxy 
model seen edge on as eSN is varied. The box is 20kpc 
across. 



length (Rd) is us ually in the range 0-25 - 0.3 3 for nor- 
mal disk galaxies (|Bizvaev fc Mitronovall2002l ). We find 
that values of c* larger than 0.05 tend to form "thick" 
disks with a Yi^/Yld ratio larger than observed, possibly 
because more stars are allowed to form in regions fur- 
ther away from the disk plane. The Rz/Rd ratio is only 
weakly dependent on eSN or /3 (runs 1 through 4, 5, 6 
& 7 in Table 1). This suggests that at our resolution 
low values of c* should be preferred to obtain thinner 
disks. Tests in S06 show c* to be weakly dependent on 
resolution when at least a few tens of thousands gas -|- 
star particles are used. However, higher c* values have 
been use d when star forming regions are individually 
resolved (iTasker fc Brvanll2005h . 

As the j3 factor is increased, more gas is affected by 
feedback and the volume ratio of the hot/cold phases 
increases (runs 6,7,8 in Table 1). In our tests we found 
that runs with the required low values of c* tend to form 
and maintain a large bubble of hot gas in the central 
region of our models, locally inhibiting star formation, 
unless /3 is also low. We thus excluded the region of 
parameter space with (3 > c* from our analysis as being 
unrealistic. 

Figure 2 shows how the gas spatial distribution of 
the DWFiso model is affected as c* and l3 are fixed to 
their desired value of 0.05 and eSN is varied from 0.1 to 
0.6. In the DWFiso runs, higher eSN leads to lower star 
formation rates and a higher cold gas turbulence (e.g 
runs 3, 4, and 8). The SFR of the dwarf galaxy model 
is most sensitive to the SN energy efficiency parameter 
eSN, with a realistic SFR obtained with values in the 
range 0.1-0.6. The images in Figure 2 show that large 
SN efficiency values create a "galactic fountain" and a 
patchy gas distribution on the disk plane resembling 
that of observed small galaxies. However, even with e = 
0.6 only a very small fraction of the gas ejected from the 
disk plane becomes unbound from the dark matter halo, 
in agreement with ob servational evidence and some the- 
oretical expectations ( Mac Low fc Ferrarall 19991 : iMartinI 
ll999l : lMaver fc Moorj|2004 lMcGaughll2005h . 

Tests on the MWiso model in S06 preferred val- 
ues of c* in the 0.01-0.1 range but did not constrain 
eSN significantly due to the relatively deep potential. 



For the MWiso model we set both c* and /3 to 0.05 
and varied eSN in the range 0.1 to 0.6 (Table 2). As in 
S06 the effect of injecting more energy into the ISM is 
much smaller than for DWFiso, due to the deeper po- 
tential well. Increasing eSN by a factor of 6 has a small 
effect on the properties of the ISM, with agas, and the 
hot/cold gas volume ratio roughly doubling. A smaller 
"galactic fountain" is created, but only a small frac- 
tion of the gas is ejected away from the disk and none 
actually leaves the boundaries of the surrounding dark 
matter halo. The vertical velocity dispersion of the stel- 
lar disk increases slightly while the SFR decreases by 
less than 20% across the whole disk. Disk galaxies over 
a range of stellar masses show a remarkably tight rela- 
tion between the gas surface density Egas and the local 
surface density o f star formation, T^sfr- Equation 4 of 
lKennicufr3 (|l99 8') (originally formulated for the average 
properties of individual galaxies) states the exact form 
of this relationship: 



(2.5±0.7)xl0~''( 



n1.4±0.15 



M^yr ^kpc 



As expected from the formulation of eq.5 



l|Kravtsov fc GnedinI [ioosh . star formation across the 
whole disk of our MWiso models follows closely the slope 
of the observed Kennicutt's law and is shown in Figure 
3. Increasing eSN lowers the SFR normalization slightly, 
again in agreement with the finding from S06. 

To summarize the set of tests described in this sub- 
section shows that our SF -I- SN feedback scheme pro- 
vides a reasonable match to some of the basic properties 
of present day, isolated galaxies: Kennicutt's law, the 
observed SFR of galaxies, the turbulence of cold gas (at 
scales of a few kpc) and the thickness ratio of stellar 
disks. This provides us with some confidence that our 
simple scheme is able to capture the important aspects 
of SF and its effects on the ISM of galaxies over a signif- 
icant range of masses. There is some degree of interplay 
between the parameters that could be completely cap- 
tured only by an even larger set of tests. However, it 
is clear that within the range of parameters explored 
in S06 and here, eSN predominantly regulates the SFR 
and the turbulence of the gas. /3 has to be relatively 
small to avoid unwanted artificial effects and c* plays a 
relevant role in creating "thin" stellar disks. Given all 
of the above considerations we have chosen to run our 
cosmological simulations keeping c* and /? fixed at 0.05 
while varying eSN in the range 0.2 - 0.6. 

3 COSMOLOGICAL RUNS 

For our cosmological runs we selected three galaxy- 
sized halos (DWFl, MWl & GALl) from 28.5 and 100 
Mpc boxes, low resolution, dark matter only simula- 
tions run in a concordance, flat, A-dominated cosmol- 
ogy: 0,0 = 0.3, A=0.7, h = 0.7, as = 0.9, shape pa- 
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Run 


c* 


/3 


eSN 


SFR 


'^gas 


^star 




Hot/Cold Gas Volume Ratio 










M0/yr 


km/sec 


km/sec 






MWisol 


0.05 


0.05 


0.1 


0.85 


12 


13 


0.04 


0.22 


MWiso2 


0.05 


0.05 


0.4 


0.8 


18 


20 


0.05 


0.26 


MWisoS 


0.05 


0.05 


0.6 


0.7 


20 


24 


0.05 


0.4 



Table 2. Summary of the properties of the isolated Milky Way runs: quantities defined as in Table 1 




0.1 1.C 10.0 



Figure 3. Star Formation rate vs local gas surface density. 
Dots from the MWisoS galaxy model run (eSN=0.6). The 
straight line is a fit to Kennicutt's Law. 

rameter F = 0.21, and = 0.039 (Perlmutter et al. 
1997, Efstathiou et al. 2002). Given the masses and 
virial radii of the selected hales the sizes of the starting 
boxes are large enough to provide realistic torques. The 
power spectra to model the initial linear density field 
were calculated using the CMBFAST code to generate 
transfer functions (Zaldarriaga & Seljak 2000). These 
three halos were resimulated at higher resolution using 
the volume renormalization technique (Katz & White 
1993). 

We resimulated DWFl, MWl, and GALl to have a 
similar dynamical range in each simulation. Each galaxy 
model was resimulated to have a similar number of par- 
ticles within the virial radius, and the softenings were 
rescaled so that the spatial resolution of each galaxy was 
a similar fraction of its virial radius. With our choices 
of particle number and softening, the smallest subha- 
los resolved have typical circular velocities of 20% of 
their host. For all particle species in the high resolution 
region, the gravitational spline softening, e{z), evolved 
comovingly from the starting redshift (z ~ 100) un- 
til z=9, and remained fixed at their final value from 



z=9 to the present. The softening values chosen are a 
good compromise between reducing two body relaxation 
and ensuring that disk scale lengths and the central 
part of dark matter halos will be spatially resolved (see 
IPiemand ct al. 2004 for a number of relevant tests). In- 
tegration parameter values were chosen as suggested in 
Moore et al. (1998) and then confirmed in Power et al. 
(2002). 

We selected a halo mass range to study halos asso- 
ciated with typical disk galaxies. The three selected ha- 
los have masses of 1.6 lO" (DWFl), 1.15 10^^ (MWl), 
and 3.1 10^^ M© (GALl) measured within their virial 
radius R„ir (the radius enclosing an overdensity of 100 
times the critical density pcrit)- GALl, the largest halo, 
is more massive than that of our Milky Way, as a re- 
cent analysis points to a Milky Way halo of about 10^^ 
Mq (jKlvpin et al.ll2002l ). This is the mass of our in- 
termediate size halo (MWl). Our least massive halo, 
DWFl corresponds to that of a typical disk field galaxy. 
The three galaxies span a range of circular velocities 
Vc (measured at Kvir) and defined as \fM{r < R)/R) 
between 70 and 185 km/sec. Note that this Vc is de- 
fined differently than the Vrot used later in the paper. 
Dark matter particle masses in the high resolution re- 
gions were 7.6 x 10^ 6.05 x lO'' and 2.3 x lO'' Mq and 
the spline softening was 0.3, 0.6 and 1 kpc for DWFl, 
MWl & GALl runs, respectively. The MWl halo was 
also run at higher resolution to test resolution effects. 
The results of the resoluti on testing are in jj 8. GA Ll is 
the same halo described in lGovernato et al] (|2004l ). The 
main halo properties are summarized in Table 3. 

The merging histories and angular momentum of 
parent dark matter halos play a major role in defining 
the final properti es of the galaxies that form inside them 
(|Cole et al.ll2000D . It is therefore important to make sure 
that our halos have merging histories and spin parame- 
ters somewhat representative of the global population. 
The three galaxies were selected with the only criteria 
of the redshift of their last major merger (zi,nm) > 2 
and with no halos of similar or larger mass within a few 
virial radii. A major merger is defined here as having 
a 3:1 mass ratio. The three halos have formation times 
(defined as the main progenitor achieving 50% of its fi- 
nal mass) in the 0.6 - 1.5 redshift range. Their spin pa- 
rameter, A (defined as JE^^'^ /GM^''^) varies from 0.01 
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Run 


Virial Mass 


Virial Radius 


Vc" 


Formation time 


A 


Last major merger 


e 


^tot at z=0 




Mo 


kpc 


km/sec 


z 




z 


kpc 


dark+gas+stars 


DWFl 


1.6 10" 


142 


70 


1.5 


0.01 


2.3 


0.3 


~ 860.000 


MWl 


1.15 10^2 


271 


134 


0.6 


0.07 


2.5 


0.6-0.3^ 


~ 700.000 - 4.000.000^ 


GALl 


3.1 1012 


380 


185 


0.96 


0.035 


2.75 


1. 


~ 480.000 



Table 3. Summary of the properties of the three cosmological halos: " Circular velocity at virial radius, number of gas and 
star particles changes slightly depending on eSN. ^ smaller/larger values are for the z=0.5 Super High Res. Run. respectively 



for DWFl to 0.05 for MWl, with th e average value for 
cosmic halos ~ 0.035 (|Gardnedl200ll ). 

7.1mm is likely a cruci al parameter in defininR th e 
properties of a galaxy disk (jSteinmetz fc Navarrol[200^ . 
If mergers efficiently destroy disks, a low zimm. would 
leave less ti me to grow a new s tellar disk from newly ac- 
creted gas ijBaugh et al.l[l99^ . We therefore stress that 
the zimm of these three halos is close to the average of 
the population of galaxy sized halos in our adopted cos- 
mology, as measured iri a recent large set of N-body sim- 
ulations (|Li et al-lbdosh . The same authors showed that 
previous estimates of the average zi^m estimates based 
on the extended Press & Schechter approach were bi- 
ased towards lower z's. DWFl has its last major merger 
at z—2.3, when a head-on encounter generates a very 
prolate halo. It only accretes relatively small halos af- 
ter that. After its last major merger at z = 2.5, MWl 
has a counter rotating minor merger at z = 2. A rela- 
tively high zimm for the MWl halo is consistent with 
the Milky Way forming several Gyrs ago (Wyse 2002). 
At z = 2.75, GALl undergoes multiple mergers and ac- 
cretes several satellites after that. A few get disrupted 
after several passages through its disk. In the future, we 
plan to simulate the formation of galaxies with a wider 
range of initial conditions. 

For each halo a sequence of runs of increasing com- 
plexity was performed. Table 4 summarizes these runs. 
For the MW halo initial conditions, we performed the 
following runs: DM only (MWldm), no gas cooling or 
SF (MWlad), no feedback and no UV (MWlgO), UV 
turned on but without stellar feedback (MWlgl), and 
three runs with various supernova efficiencies (MWlg2 
to MWlg4). In all runs with SF we kept the star forma- 
tion efficiency fixed and the mass factor /3 and varied 
only the fraction of energy from SN dumped into the 
ISM eSN from 0.2 to 0.6. Based on the results from §2 
we will mainly focus on the properties of the runs with 
a SN energy efficiency of 40% (eSN=0.4). Results from 
the cosmological runs confirm a posteriori this choice: a 
higher SN efficiency creates a DWFl galaxy with an un- 
realistic low cold gas mass fraction in the disk (< 3%). 

Our MWlg4 (eSN=0.6) run of the MWl halo was 
repeated at eight times better mass resolution and 2 
times better spatial resolution (run MWlhr in Table 4). 
The run was stopped at z=0.5 due to its computational 



run 


c* 


/3 


UV 


eSN 


DWFlgl 


0.05 




off 


off 


DWFlg2 


0.05 


0.05 


on 


0.4 


DWFlgS 


0.05 


0.05 


on 


0.6 


MWldm 


dm only 








MWlad 


no SF, no coolinj 








MWlgO 


0.05 




off 


off 


MWlgl 


0.05 


0.05 


on 


off 


MWlg2 


0.05 


0.05 


on 


0.2 


MWlgS 


0.05 


0.05 


on 


0.4 


MWlg4 


0.05 


0.05 


on 


0.6 


MWhr 


0.05 


0.05 


on 


0.6 


GALlgl 


0.05 


0.05 


on 


off 


GALlg2 


0.05 


0.05 


on 


0.2 


GALlg3 


0.05 


0.05 


on 


0.4 



Table 4. Main Parameters of Cosmological Galaxy Runs 

cost. With four million resolution elements within the 
virial radius, a star particle mass of 3 X lO'' M© and a 
force resolution of 0.3 kpc, this run has likely the highest 
resolution ever achieved for a Milky Way sized galaxy 
simulation carried to a relatively small redshift. We will 
use this simulation to test the effects of increasing res- 
olution in §8. 



4 THE EFFECT OF FEEDBACK ON THE 
PROPERTIES OF SATELLITES 

To explore the effect of UV and supernova feed- 
back on satellites, we compare the properties of satel- 
lites within the virial radius in the six MWl simulations 
described in Table 4. The MWl satellites are the most 
easily compa rable to the satellites of our own Milky Way 
(|Matedll99j '). Unless specified, similar qualitative con- 
clusions apply to DWFl and GALl. We will present 
a much more detailed analysis of the satellites in a 
subsequent paper. Sub halos were identified using SKID 
(|Governato et al1ll997l ) and only those with at least 64 
DM particles were used in the analysis. That minimum 
number of dark matter particles translates to minimum 
dark matter masses of: 0.49, 3.9, and 14.9 x 10® Mq so- 
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Figure 4. The luminosity functions of satellites within 
Hyir in all 4 MWl runs, compared with the luminosity 
functions of both the Milky Way's and of M31's satellites. 
Only satellites containing at least 64 DM particles are 
included. The 2 runs without supernova feedback produce 
far too many satellites. However, all 3 of the runs including 
supernova feedback produce a number of satellites that 
reasonably matches that of the Milky Way and M31 
satellite populations (lower and upper solid lines 
respectively), although the simulated satellites are still too 
bright. 

lar masses for the satellites of DWFl, MWl, and GALl 
respectively. The resolution limit translates into a cir- 
cular velocity limit of Vc ~ 20-30 km/sec depending on 
the individual satellite density profile. We use the ages 
and metallicities of each star particle in the satellites to 
determine their absolute V-band magnitudes assuming 
no dust reddening. Satellites have a mean (not luminos- 
ity weighted) B-V of 0.63. We verified that field dwarfs 
have slightly bluer colors. 

Figure 4 shows the satellite luminosity functions of 
the MWl cosmological runs compared with those ob- 
served for the Milky Way and M31. A comparison of 
both the gO and gl simulations (no feedback; then UV 
on, but no SN feedback) with the observed satellite lu- 
minosity functions of the Milky Way and M31 high- 
lights the missing satellite problem. On the contrary, 
runs including SN feedback result in satellite popula- 
tions similar in number, although not in luminosity, to 
the observed populations. The total number of luminous 
satellites is not strongly dependent on the strength of 
feedback implemented in the S06 algorithm. 

Although both gO and gl have a total of 26 satel- 



Figure 5. The baryon to DM ratio for all the subhalos 
within Ryir of the MWl galaxy. With no UV field or SN 
feedback the baryon fraction is constant down to our 
resolution limit of 64 DM particles. UV creates some dark 
satellites, while SN feedback removes a substantial fraction 
of baryons from halos with total masses below a few times 
lO' M0. Note that total satellites masses have been 
affected by tidal stripping. Halos with no baryons in the 
UV only run have been shifted upward for clarity. 

lites, the runs including supernova feedback all have 20 
or fewer satellites total, only some of which host stars. 
This difference in total satellite number occurs because 
simulations with more severe overcooling form satellites 
very centrally concentrated and less susce ptible to tidal 
effect s and complete disruption (see also iMaccio et al.l 
I2OO5I ). Similarly, runs with cooling but no SN feedback 
(as MWlgO and MWlgl) contain ~2x the number of 
satellites as the DM only run or the run with no gas 
cooling. 

Figure 5 shows that some of the satellites in the 
runs including UV are completely devoid of stars. There 
are no satellites within Kvir in any of the simulations 
that retain gas but don't contain any stars. Unlike UV 
feedback, supernova feedback can also expel gas from 
the satellites once they start forming stars, reducing 
the overall baryon fraction of satellites containing stars. 
These reduced baryon fraction satellites are seen in Fig- 
ure 5. As expected, the satellites of the run with stronger 
supernova feedback have a smaller baryon fraction than 
the run with weaker supernova feedback. The luminos- 
ity functions of satellites in the MWl (and of GALl) 
simulations show that the introduction of SN feedback 
in particular is essential to reproduce the observed num- 
ber of satellites in Milky Way type galaxies. The number 
of satellites containing stars appear is fairly robust as 
resolution is increased (see Fig. 20). At a higher reso- 
lution while tens of less massive dark satellites are re- 
solved, only a few more are able to make stars. This is 
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Figure 6. Left panel: Brightness maps of the edge-on disks of (from top to bottom) GALl, MWl and DWFl at z=0. Each 
star particle in the simulations has been weighted by its age-dependent bolometric luminosity. Right Panel: the face-on surface 
density of the gas for DWFl at z =0. Each frame is 30 kpc across 



because the mass threshold belo w which the UV fi eld 
makes satellites completely dark (jOuinn et al.l[l996l ) is 
resolved in our runs (Fig. 5). All but one satellite of 
galaxy DWFl (the run with the best mass resolution 
in our set) are completely dark when SN feedback is 
on. Note that the list of known MW satellites at very 
low luminosities and surface br i ghtne sses might still be 
incomplete (see lWillman et aP ( 20041 ) for a discussion). 

5 PROPERTIES OF THE DISKS 

The images in Figure[6]show that all three galaxies have 
formed a significant stellar disk by redshift z = 0. In this 
section, we discuss in detail the spatial and kinematic 
properties of the disks, investigate the effects of feedback 
on disk properties, and demonstrate the formation of the 
disks during major, gas-rich mergers. 

5. 1 Density Profiles of the Disks 

The top panel of Figure [7] shows the z = stellar sur- 
face density profiles of the disks as a function of ra- 
dius for our fiducial runs (DWFlg2, MWlgS, GALlg2 - 
eSN=0.4). Hereafter we will simply refer to these specific 
runs as DWFl, MWl, and GALl as we report and dis- 
cuss our results. We include all stars within 4 disk scale 
heights, Rz. The stellar disk component of all galaxies is 
well-fit by an exponential distribution between one and 
three disk scale lengths. The profiles in Figure [7] illus- 
trate that all galaxies also show evidence for a central, 
steeper component that extends to 1 or 2 kpc. 

Disk scale lengths, Rd, were measured by fitting an 
exponential profile to the disk component, but exclud- 
ing the steep central component and stopping where the 
surface density shows a break. To approximate the effect 
of measuring the disk scale length in different optical 



bands, we measured the scale lengths of stars in differ- 
ent age ranges: < 2 Gyrs, < 4 Gyrs and < 11 Gyrs, to 
approximate the B, R and K optical bands, respectively 
l|Martin et al.l |2005| ) . These scale lengths are summa- 
rized in Table 5. The bottom panel of Figure [7] shows 
the stellar density profiles of MWl in each of these 
three stellar age ranges. Beyond 2.5-3 scale lengths 
(shown by the outer vertical line), the disk surface den- 
sity declines rapidly for the youngest stars. These disk 

J irofiles are qu ite comm on in obs erved spiral galaxies 
Pohlen fc TruiillQ .200d: iFlorido et al. 2001). In the re- 
mainder of this section we will refer mostly to stars 
younger than 4 Gyr (R band) to illustrate our results, 
unless otherwise specified. 

The exponential disk component forms from 
the inside out as ob served in many disk galaxies 
l|Rvder fc Dopitalll994l ). However, gas inflow due to the 
late formation of a bar or an oval-like structure in 
the central few kpc, which is responsible for the ob- 
served central steepening of the profile, can trigger en- 
hanced star formation in the inner disk, which then 
forms partially outside-in. This mode of disk formation 
is significant in ru n GAL l and was already reported 
in iGovernato et ahl (|2004h . At z = 0, star formation in 
GALl is primarily concentrated in the inner parts of 
the disk. The scale length of stars younger than 2 Gyr 
in GALl is only 1.1 kpc, half that of the older com- 
ponent. This is typ ical of massive SO/early type disk 
galaxies (| Pohlen et al. 2004). 

All three galaxies have extended disks of gas colder 
than a few 10'* K with exponential scale lengths larger 
than even the y o unges t stellar disk component, as in 
iRvder fc Dopital (|l994h . The right panel of Figure [6] 
shows that at our sub-kpc resolution and with our feed- 
back scheme we are able to qualitatively capture the 
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Figure 7. The face-on stellar density profile of the three 
simulated galaxies. Only stars within 4 Rz from the disk 
plane are considered. The upper panel shows the 
contribution from all stars with a cosmic age between and 
10 Gyrs, roughly co rresponding to the K band 
jMartin et al.ll2005l) . The lower panel shows the projected 
stellar density profile of the MWl galaxy for stars in the 
0-3 Gyrs (solid), 0-4 Gyrs (dotted) and 0-10 Gyrs (dashed) 
age range. Dots show the stellar profile of the MWlhr at z 
= 1.7. The two straight lines show the slope of the younger 
stellar component inside and outside the break radius. 

complex nature of the ISM. The gas distribution in the 
disk plane of our smallest galaxy (DWFl) and in the 
high resolution MW (MWhr) shows a complex struc- 
ture of spiral arms linked by filaments and surrounded 
by a hotter medium. 

The stellar profile of the final disks in the more 
massive galaxies (MWl and GALl) is only weakly de- 
pendent on the strength of supernovae feedback as we 
allow the energy efficiency to vary only by a factor of 
thr ee, being constrained b y the early tests we carried in 
§2. lOkamoto et ahl (|2005h showed how the density dis- 
tribution might depend substantiaUy on the feedback 
adopted. It is likely that the energy injection per unit 
mass in the feedback scheme adopted by Okamoto et 
al. is higher than in our simulations. Instead, the disk 
stellar profile of DWFl, our less massive galaxy, is more 
affected when varying the amount of SN feedback in the 
same range (Fig. 8). When e increases to 0.6 the stellar 
density profile of younger stars in DWFl is best de- 
scribed by a single, exponential profile with a shorter 
scale length and truncated at about 2 R^. Similarly to 
the tests performed in §2 the cold gas is much more 
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Figure 8. The face-on stellar density profile of the DWFl 
galaxy as a function of feedback strength. Only stars within 
4 from the disk plane and in the 0-4 Gyrs age range are 
plotted. 

turbulent, making star formation inefficient at low den- 
sities. When feedback is turned off cold gas settles in a 
smaller, more compact disk that originates a denser cen- 
tral core. As the disk surface density increses the disk 
becomes bar unstable redistributing the stellar material 
and originating the "flat" stellar profile at ~ 2kpc that 
is usually associated wi th strong dynamical instabilities 
(|Debattista et al]|2005l ). In all three galaxies runs with- 
out feedback form "hotter" stellar disks with a higher 
scale lenght. 

5. 2 Kinematics of the disks 

We first studied the dynamical properties of the sim- 
ulated objects and verified that the disks are indeed 
supported by rotation. For every star particle within 
the "disk" we computed Jz, the component of angular 
momentum parallel to the total angular momentum of 
all disk stars. We defined the disk stars as those within 
4 (R band) scale lengths and heights centered on each 
galaxy. We compared J along the z axis to the angu- 
lar momentum of the co-rotating circular orbit with 
similar orbital energy, Jc(E). Stars in a disk compo- 
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Figure 9. Mass weighted probability distributions of the 
orbital circularity J/Jc for disk stars in DWFl (solid), 
MWl (dotted) and GALl (dashed). The right panel (note 
the change of scale) shows only star particles formed in the 
last 3 Gyrs. 



nent completely supported by circular motions should 
by definition have J/Jc = 1. Figure |5] shows the mass 
weighted distribution for all the disk star particles (left 
panel) and for disk star particles formed in the last 3 
Gyrs (right panel). The histogram distributions show 
clearly that in all three galaxies the disk component is 
mostly supported by rotation. Our DWFl and MWl 
galaxy models have a higher fraction of stars on circu- 
lar orbits when c ompar ed wi th the distributi o ns sho wed 
by lAbadi et al.1 l|2003l ) and lOkamoto etHI (|2005l ). In 
particular the young stellar component is heavily dom- 
inated by stars on almost circular orbits: 80% of young 
stars in the DWFl run have J/Jc > 0.8. 62% of all stars 
in the MWl galaxy have J/Jc > 0.8. As a compari- 
son most thin disk star s in the MW have J/Jc > 0.7 
l|Nordstrom et ai]|2004l ). A central hump in the J/Jc 
distribution, that has often been interpreted as the sign 
of a massive halo/bulge component is evident in the 
DWFl run. We verified that the low angular momen- 
tum stars are also the oldest, and were formed during 
DWFl's last major merger. 

The amount of stellar rotation versus velocity dis- 
persion, w/cr, is another good measure of how realis- 
tic are the disks that we obtained: Vrot/cr ~ 4 for the 
MWl disk stars younger than 4 Gyrs, where a is the 
rms average between the tangential and radial velocity 
dispersions. Typical observed values range from 2 to 5 
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Figure 10. The V/ct ratio for the young stellar component 
of DWFl as feedback is increased: Long dashed: no 
feedback or UV, dot dashed: eSN = 0.4, solid: eSN = 0.6 



for the disks of normal galaxies ijPizzella et al.ll2004l ). 
Fig. 10 shows Vrot/o" for the young stellar component of 
DWFl. The dynamically coldest stellar disk is obtained 
with eSN = 0.4. Stronger feedback makes the ISM more 
turbulent, causing stars to form with a higher velocity 
dispersion. For all three galaxies runs without SN feed- 
back generate stellar disks with lower N rot jo ratios, due 
to stronger loss of angular momentum driven by bar in- 
stabilities or because the velocity dispersion is increased 
by the overall stronger heating by satellites (that sur- 
vive disruption much closer to the galaxy center) . These 
statements hold true irrespective of the selected age in- 
terval for the stars. 

The azimuthally averaged rotation curves of galax- 
ies are shown in Figure 1111 Both the rotation curve of 
the cold gas component and of the disk stars younger 
than 4 Gyrs are plotted. For both the MWl and DWFl 
runs intermediate age stars and cold gas have very sim- 
ilar profiles, confirming that younger stars are in orbits 
close to circular. In the GALl run (our most massive 
galaxy) , the stellar and gaseous disks are not completely 
coplanar. This offset creates a difference in the central 
part of the velocity profile of the two components. The 
outer part of GALl's stellar disk is also dynamically 
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Figure 11. Azimuthally averaged rotation velocity of cold 
gas (dot — dashed line) and stars ( age < 4 Gyrs). Prom top 
to bottom: GAL1,MW1,DWF1. The vertical lines show 
2.2 Rjj for each galaxy. 

hot, which makes Vrot of the stellar component decline 
very rapidly at outer radii. 

While all our galaxies show long lasting bar in- 
stabihties, the introduction of SN feedback increases 
the stabihty of the disk and thus reduces the non- 
axisymmetry of the stellar disk. Stronger feedback gen- 
erates stronger disk axisymmetry. In all simulations 
with feedback, only an oval distortion with scale smaller 
than the disk scale length is present as opposed to a clear 
bar with scale length ~ disk scale length in models with- 
out SN feedback. The small oval distortion is likely an 
unresolved bar, because our high resolution control run 
(MWlhr) does display a resolved bar at its center. Disks 
with feedback are more stable to bar formation because 
their stellar disk builds up more slowly over time. 

We can understand the stabilizing effect of feedback 
by comparing the mass in the halo to that in the disk us- 
ing th£_gaxameter£ X Rd) / GMd, introduced 
by lEfstathiou et all (|l982l ) to measure the susceptibil- 
ity of disks to bar formation. Vmax is the peak circular 
velocity and Md and Rn are, respectively, the mass and 
exp one ntial scale length of the disk. Numerical experi - 
ments (lEfstathiou et al.lll982l : [Maver fc Wadslevll2004 ) 
show that e > a, with 0.95 < a < 1.1, is required for 
stability against bar formation of a disk of gas and stars 
(the higher the gas fraction, the higher q). Measuring e 



in a cylinder within 0.5 kpc of the disk plane, we find 
that in the runs without feedback e is already < 0.9 at 
2 > 1. In the runs with feedback, e does not approach 
unity until near 2; = 0. In runs without feedback the 
disk is more massive and denser at earlier epochs, and 
thus more bar-unstable. Without feedback, all the gas 
that is accreted rapidly makes disk stars, while with 
feedback a significant fraction of the baryons remain 
hot and diffuse and do not form stars even when they 
reach the disk. A visual comparison with real rotation 
curves shows that our simulated galaxies produce rota- 
tion curves that raise rapidly and then decline too fast 
at R > 4 Rfj. Only DWFl has a rotation curve close to 
being "fiat" as observ ed for most normal spiral galaxies 
(|Catinella et al.ll2006l ). This is likely due to a combina- 
tion of a still too massive central component and the 
outer part of disks being affected by two body heating 
from halo particles. 

5.3 Gas Rich Mergers and the Formation of Disks at 
High Redshift 

With our relatively high resolution we are able to fol- 
low the formation of disks early on in our simulations. 
Perhaps surprisingly, the disk of the MWl galaxy starts 
forming immediately after the last major merger event 
confi rming early results by (|Springel fc Hernauistll2005h 
and (|Elmegreen et aLll2005l ) that gas rich mergers with 
substantial feedbac k are conducive to th e observed early 
formation of disks (jTruiillo et al.l [20051 ). In our simula- 
tions feedback stops a substantial fraction of gas from 
turning into stars during the major merger. 40% of the 
mass of the disk of MWl is gas at z = 2. With lots of 
orbital angular momentum available from the merger, 
gas cools and rapidly settles into an extended exponen- 
tial disk, especially in runs with stronger feedback and 
higher resolution. Even for DWFl and GALl, substan- 
tial stellar disks form immediately after zimm- Without 
feedback a much more compact stellar component forms 
and becomes rapidly non-axisymmetric, as described in 
§5.2. At z>l stellar disks are 20-30% larger when feed- 
back is on, although by 2 = their sizes are virtu- 
ally identical. V/cr measured in the disk of the high 
resolution MWlhr run is ~ 1.6 a.t z — 1.5 for stars 
that form shortly after zimm- This stellar component is 
heated by a combination of the rapidly changing cen- 
tral potential and by a couple of minor satellite ac- 
cretions onto the disk that occur before z =1.5. This 
early disk could be associated with th e thick disk com- 
ponents of present day spiral galaxies dSeth et al.l l2005l : 
lYoachim fc Dalcantonll2005l : iBrook et al.ll2004h . 

5.4 Bulge & Disk Decomposition 

To separate the bulge, disk and halo components we 
assigned to the "bulge" all stars in the age range that 



15 



© 0000 RAS, MNRAS 000, 000-000 



14 Governato et al. 



Run 




Rdfl 


RdK 


Rzfl 








Itot 




SFRl M0/yr 






DWFl 


2.02 


1.95 


1.2 


0.3 


7. 


0.4 


-19.8 


-21.2 


-22.6 


0.6 


160 


1.38 X 10^'' 


MWl 


4.5 


3.8 


2.1 


1. 


10 


1. 


-21.1 


-22.6 


-24 


2.2 


205 


4.6 X 10i° 


GALl 


1.1 


1.2 


2.15 


2.3 


6. 


1.3 


-22.0 


-23.5 


-24.9 


7.1 


180 


1.36 X 10" 



Table 5. Summary of the properties of the disks of the galaxies formed in the cosmological runs (only runs with eSN = 0.4 
are shown). ^ at z=0, ^ at 2.2 R^j. 



would give a Nrotja < 1 and within four disk vertical 
scale heights from the disk plane. This selects an old, 
slowly rotating stellar component. Disk stars were those 
within four disk scale heights and scale lengths that were 
not assigned to the bulge. Stars not belonging to the 
bulge, disk or individual satellites were assigned to the 
halo. The kinematically defined B/D mass ratio of the 
fiducial MWl is = 0.3 or smaller. This result does not 
change for other reasonable decomposition techniques 
l|Governato et all [20041 . We only used dynamical def- 
initions of the galaxy components. A photometric de- 
composition was not attempted but given that stars in 
the disks are younger and then brighter it is likely that 
a luminosity weighted estimate would give a lower B/D 
ratio. We speculate that the low B/D ratio of the MWl 
is due to the high angular momentum acquired by cold 
gas during the last major merger event. The DWFl and 
GALl simulations both formed a more massive "bulge". 
This result is likely a consequence of different physical 
processes: (i) The DWFl galaxy has very low angular 
momentum (the last major merger at z = 2.3 is almost 
head-on); (ii) GALl forms stars in its central region 
until the present time, as gas keeps cooling from the 
surrounding halo. Its bulge has fairly blue colors and is 
relatively more massive. The relevant scale lengths are 
< 1 kpc, so we decided to postpone a more detailed anal- 
ysis of the bulge components of our simulated galaxies 
when higher resolution simulations will be available. 



MATCHING THE TULLY 
LATION 



FISHER RE- 



The TuUy-Fisher (TF) relation links the characteristic 
rotatio n velocity of a galaxy w ith its total absolute mag- 
nitude jGiovanelli et al. 1997 ). Similarly, the "baryonic 
TF" relation (|McGaugh et al] I2OO0I : iMcGaughl |2005| ) 
links characteristic rotation velocity with total disk 
mass to account for the fact that less massive galax- 
ies are more gas rich and thus that their stars only 
account for a small fraction of their disk total mass. 
Both relations imply that the total disk mass and its 
radial distribution are closely connected with the to- 
tal mass of a halo. Early simulatio ns matched the slope 
of the TF relation, but reported l|Eke et al-lboOll ) dif- 
ficulties in matching its normalization. The relatively 
weak feedback adopted in that work might have caused 



galaxies to cease star formation relatively early, leav- 
ing the disks too red and then causing the TF relation 
offset. Alternatively, the central rotational velocities of 
the simulated galaxies were too high due to an excess 
of matter at the center of galaxies possibly due to ex- 
cessive baryonic cooling and the subsequent adiabatic 
contraction of the dark m atter component . Even high 
resolution simulati ons like iGovernato et ah I l|2004 ) and 
lAbadi et al.1 (|2003h have shown a consistent offset from 
the TF relation of late type spirals, with galaxies be- 
ing too centrally concentrated compared to real ones. 
Compared to some previous works our simulations have 
sufficient resolution to allow us a much direct compar- 
ison with observational data as Vrot is now measured 
at only a few kpc from the center. This poses a much 
stronger constraint on simulations as they have to re- 
produce the stellar mass distribution and kinematics at 
relatively small scales. 

When feedback is included, simulations successfully 
match both the TF and baryonic TF relations, this re- 
sult being possibly the most important conclusion of the 
work presented here: it shows that our simulations pro- 
duce a realistic distribution of stars, baryons and DM 
within the central few kpc. To compare our simulations 
with the observed TF relations, we used the rotation 
curves plotted in Figure 9 (and described in §5.2) to 
determine Vrot measured at 3.5 disk scale lengths for 
DWFl, MWl, and GALl as in the GiovaneUi sample. 
We then obtained the global magnitudes for each galaxy 
model by coupling the star format ion histories (SF Hs) 
of our simulations with GRASIL (|Silva et al.lll99^ . a 
code to compute the spectral evolution of stellar sys- 
tems taking into account dust reprocessing. We con- 
sidered each star particle as a single stellar population 
with the age, metallicity, and mass assigned to it by the 
simulation. Stars with metallicities less than Z — .0005 
were assigned that minimum value. GRASIL models the 
distribution of molecular clouds and dust based on the 
structural parameters of a galaxy disk and bulge. The 
quoted magnitudes in Table 5 are the average of differ- 
ent viewing angles. 

Figures [12] and [13] show our galaxy models over- 
plotted on TF and baryonic TF da ta from Gi o vaneU i 
(2005, private communication) and iMcGaughl (|2005l '). 
respectively. The opposite residuals for DWFl and 
MWl from the average TF relation (continuous line in 
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Figure 12. The TuUy Fisher relation using the data 
compilation from GiovaneUi 2005 (private communication) 
and a fit to lGiovanelli et alj lll997t) . Solid triangle: DWFl, 
Solid Square: MWl, Solid Dot: GALl. Bigger dots shows 
V, ot measured at S-SR^j. Smaller dots shows the efli'ect of 
measuring Vrot at 2.2Rrf. The small open dot uses Vrot 
measured from GALl cold gas component 



Figure 13. The Bary onic TuUy Fisher relation. Squares 
show data points from l lMcGaugh 2005). Vrot is measured 
at 3.5Rd (as in the observational sample) for each of the 
three galaxies. The two smaller points for the GALl run 
show the effect of measuring Vrot at 2.2 R^ using stars and 
the cold gas component, as shown in Figure [12] 



Fig. 11) correlate well with their halo spins: lower than 
average for DWFl and higher for MWl. As GALl dy- 
namical tracers are slightly decoupled (see §5.2), its lo- 
cation on the TF relations somewhat depends on the 
tracers chosen (old or young stars for the disk scale 
length, stars or cold gas as tracer of the velocity field). 
As GALl's stellar mass is larger than lO'^^ M0 it is a 
likely candidate for an early spiral or even an SO. We 
speculate that another type of feedba ck, possibly from 
an a c tive QSO at the b ulge center (|Di Matteo et all 
I2OO5I : iHopkins et aLllioOsI ). would decrease star forma- 
tion in the center and lower the central gas density, re- 
ducing the relative weight of the central baryonic con- 
centration. 

Improvements in matching the TuUy Fisher and 
in conserving the angular momentum content of 
disks have been consistently reported in the lit- 
erature as the robustness of simulations and the 
detail of physical modelling has continued to im- 
prove (ISommer-Larsen et aLfbOOSl: lAbadi et all l2003l : 



iGovernato et all |2004 iRobertson et all 1200^ 1 As we 
show in §8 our simulations have a softening of the order 
of only 20% of the disk scale lengths and give robust 
and likely converging results at scales corresponding to 
the typical scale length of stellar disks. We verified that 
our match with the TF relations is not due to Vrot << 



Vc, i.e. the disks are not dynamically "hot" within a 
few disk scale lengths and gas and young stars trace 
very closely the underlying potential. We verified that 
runs with no feedback produce galaxies slightly redder, 
but only in the B band and only for DWFl where the 
effect is 0.3 mag with eSN = 0.4 and 0.5 mag for eSN 
= 0.6. For more massive galaxies and/or redder bands 
the effect is only 0.1 - 0.05 mag as the effect of feedback 
is smaller and star formation happen mostly at higher 
redshifts. This small color shift is not sufficient to ex- 
plain the shift from the I band TF measured in some 
previous simulations. 

7 FEEDBACK AND THE STAR FORMA- 
TION HISTORIES OF GALAXIES 

Not only does feedback affect the structure and stability 
of disks, as discussed in §5.1 and §5.2 respectively, but 
we also expect it to shape their star formation histories. 
The isolated test simulations of §2 hinted at this effect, 
and in this Section we explore the effect of feedback 
on the star formation histories of DWFl, MWl, and 
GALl. The SFHs of bulge and disk stars in these three 
runs are shown in Figure 14. The plot shows that the 
peak of star formation is delayed in less massive halos, 
as expected if feedback has a more pronounced effect in 
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Figure 14. Star formation histories of DWFl, MWl, and 
GALl, including all stars within 4 disk radial and vertical 
scale heights. eSN=0.4 and UV field on for all runs. Solid 
black: MWl, long dashed: GALl, dotted: DWFl. 

galax ies with shallower potential wells (|Maller fc Dekell 
l2002h . To understand this result, consider the natural 
property of CDM models that star formation at high 
redshift happens in a number of progenitors that even- 
tually coalesce to form a present day galaxy. We used 
the DWFl simulation to verify that increasing feedback 
limits the efficiency of star formation in progenitor halos 
with Vc < 30-50 km/sec and delays it until small sub- 
halos merge i nto halos with a suffic iently deep potential 
weU (see also iNeistein et al.l (|2006l ')'). 

To illustrate the crucial effect of feedback on the 
SFH of DWFl, Figure [13 compares its SFH in the ab- 
sence of feedback and with eSN = 0.6, the strongest 
feedback adopted in our study. Early star formation is 
significantly reduced in the case of feedback, and we no 
longer see the initial strong starbursts evident in the run 
with no feedback. Lowering the efficiency of star forma- 
tion during mergers preserves a large quantity of cold 
gas that rapidly settled in rotationally supported disks. 
The run with feedback has a present day SFR almost 
ten times higher than that without. The SFR in the run 
with feedback is far too low to match observations. 

The z = star formation rates of DWFl, as well 
as of the other two galaxies, are quite similar to those 
observed in galaxies with the same stellar mass. Fig- 
ure [16] shows the z = star formation rates of our 
simulated galaxies overplotted with those measured for 
SDSS galaxies in a wide mass range (jBrinchmann et ahl 
|2004| ). Two low mass field galaxies outside the virial 
radius of galaxy MWl (but within the high resolution 
region) have also been included to show a fair agree- 



Figure 15. Galaxy DWFl: SFH including all stars within 
4 Rz and R^j from the disk plane for two runs. Solid line: 
eSN=0.6, long dashed: no feedback, no UV. The addition of 
feedback smooths out the SF peaks otherwise present at 
high redshift and during the last major merger event at 
z=2.3. Feedback delays the conversion of gas into stars until 
gas accumulates and cools in the potential well of the main 
progenitor. 




-3 



8 9 10 11 

Log M,/Mg 

Figure 16. Present day star formation rate vs galaxy 
stellar mass for both simulated and observed galaxies. 
Triangle: DWFl, Square: MWl. Solid Circle: GALl. The 
two smaller dots are field galaxies in the high resolution 
region of our simulations. The continu ous line is the avera ge 
SFR f or the SDSS sample described in iBrinchmann et al.l 
l|2004h . The central dashed line is the median value. The 
upper and lower dashed lines include data three sigma from 
the mean of the observational sample 
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ment between our theoretical model and observational 
data over more than two orders of magnitude in stellar 
mass. 

Our galaxy models run with SN feedback also 
rep roduce the observed "downsizing" of fie ld galax- 
ies (jCowie et al.lll996l : iMacArthur et al.ll2004l '): galaxies 
with a less massive stellar component are also younger. 
Figure [17] shows that the inclusion of feedback and its 
differential effects at different galaxy masses naturally 
reproduces the observed trend: more massive galaxies 
also have older stellar populations. This trend disap- 
pears when runs with no feedback are considered (open 
circles) and average stellar ages correspond to the time 
of collapse of the first protogalaxies. As SN feedback re- 
moves baryons and delays star formation in smaller ha- 
los, it has a stronger effect in less massive galaxies that 
formed from even less massive progenitors at high red- 
shift. Gas must collect inside a sufficiently deep poten- 
tial well before it is able to cool efficiently. More massive 
galaxies are formed from the merging of more massive 
progenitors, that were able to form stars already at high 
redshift. While CDM forms structure from the bottom 
up, the introduction of feedback reverses the trend for 

stellar populations. 

Recently, iBell fc de Jond (|200ll ) and 
IConselice et al.l (|2005l ) used observations of spiral 
galaxies to measure the stellar mass associated 
with galaxies of different rotational velocities and 
halo mass. A comparis on of our simulations with 



data in jConselice et al.1 Hooi) and the best fit 



to 

iBeU fc de J ong (200^) is plotted in Fig.18 and shows 
that our models form the right amount of stars in their 
bulge and disk components. We also verified that MWl 
contains 4 x 10^" M0 of dark matter within 8.5 kpc, 
again in good agreement with estimates for the Milky 
Way (Eke et al. 2001). These results confirm that our 
modelling of SF and feedback created galaxies with the 
right amount of stars for their halo mass. 

8 RESOLUTION TESTS: THE CENTRAL 
MASS PROFILE OF GALAXIES AND THE 
LUMINOSITY FUNCTION OF SATEL- 
LITES. 

In this section, we show that our results are robust to 
resolution effects. To do this, we have performed one ad- 
ditional simulation (MWlhr) to test if our results have 
converged as resolution is further increased. A more 
complete study on how resolution affects the formation 
of disks and the general properties of galaxies is highly 
needed but beyond the scope of this paper. MWlhr uses 
the same initial conditions as the MWl runs, but uses 
eight times as many particles within the high resolution 
region. eSN = 0.6 (as for run MWlg4) was used. Force 
resolution is improved by a factor of two to 0.3 kpc. Ad- 
ditional waves in the initial perturbations spectrum are 
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Figure 17. Mass weighted average age of stellar 
populatio ns vs galaxy stellar mass . Small dots with error 
bars: from IMacArthur et al. Solid large dots; 

simulations with eSN = 0.4. Open dots: simulations with no 
feedback or UV. Vrot is measured at 3.5 disk scale lengths. 
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Figure 18. Total stellar mass vs. rotational ve locity. Solid 
small dots are data from lConselice et al. Large dots: 

simulations measured at 3.5 R^j. Small points: simulations 
measured at 2.2 R^. T he straight line is a fit from 
iBell fc de Jon3 1I2OOII) . 
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consistently added. With a dark matter particle mass of 
only 7.6 x 10^ Mq and star particle mass of 3 x 10'' 
Mq this simulation has one of the highest resolutions 
for a Milky Way sized galaxy in a cosmological context 
to date. Due to its high computational cost we carried 
this simulation only to z = 0.5, using several tens of 
thousands of CPU hours. Figure [19] compares the circu- 
lar velocity Vc of MWlg4 and MWlhr (again defined as 
M{r < R)/R) for the DM, baryonic and stellar com- 
ponents, while Figure[20]compares the V-band satellites 
luminosity function (again at z=0.5). 

Figure [T^ shows that at the resolution of our stan- 
dard runs the total mass profile and the quantity of 
stars have converged at radii corresponding to a cou- 
ple of disk scale lengths. Vrot measured at 3.5Rd (and 
used in Figures 12, 13, 17 and 18) is thus robust to res- 
olution effects, at the resolution of our standard runs. 
Test runs in a cosmological context resolution (Gover- 
nato et al in prep.) show a systematic increase in Vc 
and a decrease of stellar R_d as resolution becomes in- 
creasingly worse. This result likely explains why our 
simulations match the normalization of the TF rela- 
tion. However at smaller radii (only 3-4 times the grav- 
itational softening) our standard run MWlg4 shows a 
much larger central mass than MWlhr. This unphysi- 
cal mass concentration significantly raises the Vc such 
that it peaks at over 300km/sec. Because the standard 
and high resolution runs adopted identical feedback and 
SF parameters, this result shows that resolution plays 
an important role in the central region of galaxies. In 
DM-only simulations, low resolution artificially lowers 
central DM densities. In contrast with DM only simula- 
tions, low resolution and lack of feedback in SPH sim- 
ulations causes baryons to become more centrally con- 
centrated due to partially artificial angular momentum 
loss. As resolution is increased, reduced artificial angular 
momentum loss by both the stellar and gaseous compo- 
nent allow SN feedback to be more effective and to re- 
duce the central baryon concentration. Worse resolution 
and possibly weaker feedback would have overestimated 
the amount of baryons within 2.2 R^. As the force reso- 
lution is increased the disk bar instability in the disk of 
MWlhr starts earlier (z ~ 2) and is more pronounced 
than in the MWlg4 simul ation, confirming studies by 
iKaufmann fc Maverl (|200d ). At this high resolution the 
fraction of cold baryons is 40% higher, and the mass in 
stars is only 5% lower. We also begin to resolve the ther - 
mal instabilities predicted by iMaller fc Bullockl (|2004l ') 
and recently detected i n simulations of isolated galaxies 
IIKaufmann et al ] |2004h . It is encouraging that the satel- 
lite V-band luminosity function (Fig. 20) remains sub- 
stantially unaltered as resolution is increased, although 
a few more satellites as faint as Mv = —8 are now re- 
solved. This supports results from S06 that our adopted 
SF-I-SN algorithm is relatively insensitive to resolution. 
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Figure 19. Resolution tests: The circular velocity Vc= 
^/M{< R)/R for run MWlhr (solid) and MWlg4 (dots). 
The mass distribution converges at radii larger than r > 2.2 
Rd (vertical line). The high resolution simulation contains 
significantly less baryons within a few times e 
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Figure 20. Resolution tests: The V - band luminosity 
function of the satellites system of MWlg4 (dashed) and its 
high resolution version (dotted) at z = 0.5 compared with 
the Milky Way and Andromeda (solid lines). 
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9 CONCLUSIONS 

We have presented a new set of simulations of the for- 
mation of disk dominated galaxies in a full ACDM cos- 
mological context carried to the present time. These 
simulations included a simple but physically motivated 
recipe for star formation. We used isolated galaxy mod- 
els to calibrate the free parameters of our star forma- 
tion/feedback scheme to reproduce several important 
properties of the stellar and ISM components of local 
disk galaxies over a range of masses. We then used this 
SF/feedback algorithm in cosmological simulations with 
no further adjustment. 

We simulated three different halos in a cosmolog- 
ical context, with masses typical of those associated 
with disk galaxies. These cosmological simulations have 
a very high number of resolution elements within the 
virial radius of the simulated galaxies (~ 1.5 10^ for 
each of the DM and gas components and a few hundred 
thousand star particles), allowing us to resolve individ- 
ual satellites down to Vc ~ 20% of the parent galaxy. A 
rotation dominated stellar disk galaxy with a significant 
exponential component naturally forms within each of 
these halos. Our MWl simulation provides a good ex- 
ample of how the stellar disk can form shortly after a 
major merger event, due to the large amount of orbital 
angular momentum transferred to the gas component 
during the merger event. Our three simulated galaxies 
(DWFl, MWl, GALl) share the size, stellar mass, col- 
ors and current SFR typical of a low mass galaxy, of a 
Milky Way like galaxy and a massive early type spiral, 
respectively. 

One of our most significant results is that the simu- 
lated galaxies fall on the TF and baryonic TF relations. 
This match is possible in part due to the high resolution 
of the simulations: Resolution is crucial to avoid arti- 
ficially losing th e angular momentum gained this way 
(|Governato et al. 2004: Ka ufmann fc Maver 2006). At 
our resolution, rotation curves have been resolved down 
to a fraction of the stellar disk scale lengths and (mea- 
sured at 3.5 Rd) range from 135 to 270 km/sec in am- 
plitude. The MW sized galaxy contains an amount of 
DM within the solar radius well in agreement with ob- 
servational constraints. The simulated galaxies also re- 
produced the observed 'anti-hierarchical' trend of stel- 
lar populations: smaller galaxies have younger stellar 
populations. As feedback delays star formation in pro- 
togalaxies a larger reservoir of gas remains available to 
form disks until the bulk of the mass of galaxy is assem- 
bled, typically between redshift 3 and 1. Feedback and 
the cosmic UV field also drastically reduce the number 
of galaxy satellites containing a significant stellar pop- 
ulation. 

We have complemented our standard runs with a 
high resolution test run of our MW sized halo. This run 
has and unprecedented number of particles: several mil- 



lion within Kvir- We used it to test the convergence of re- 
sults in our standard runs. At our highest resolution the 
main galaxy is significantly less centrally concentrated, 
reducing its bulge mass and its peak velocity consider- 
ably. While stellar bulges are likely still unresolved in 
our standard runs, our tests and estimates from previ- 
ous works suggest that the mass distribution and the 
rotation velocity of the cold gas and young stellar com- 
ponents have converged at the scales of a few kpc. Such 
a high resolution coupled with a physically motivated 
description of feedback explains why our models do in- 
deed fall onto the TF relations. 

These simulations and those published in the re- 
cent literature show that some significant progress has 
been done toward understanding galaxy formation in 
the context of a ACDM concordance cosmology. How- 
ever, areas remain were progress needs to be made: (i) A 
larger range of initial conditions needs to be explored to 
understand the role played by halo spin and merging his- 
tory on the morphology and colors of galaxies; (ii) Our 
resolution test shows that the inner mass distribution of 
the galaxies, including the bulge, is still poorly resolved 
by our standard runs and that possibly several million 
resolution elements per galaxy are needed to approach 
convergence; (iii) Satellites are still too bright compared 
to those orbiting our Milky Way. While cosmic variance 
likely played a role (the largest satellites accreted by the 
MWl system is more massive than the LMC), one ob- 
vious possibility to create fainter satellites would be to 
increase the efficiency of SN feedback, but that needs to 
be carefully calibrated with the other constraints. Two 
intriguing possibilities that we plan to explore in future 
works are the increas ed feedback and UV field d ue to 
massive POPIII stars (jChoudhurv fc Ferraral [20051 ) and 
the local proximity effect that star formation from the 
main galaxy might have on the nearby dwarfs (Gnedi^ 
I2OO0I ) . Overall, the possibility of obtaining a good match 
to the observed satellite LF seems within reach of future 
simulations. 
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